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Abstract 

Complete two-loop calculation of a dimensionally regularized Wilson loop with light- 
like segments is performed in the light-like axial gauge with the Mandelstam-Leibbrandt 
prescription for the gluon propagator. We find an expression which exactly coincides with 
the one previously obtained for the same Wilson loop in covariant Feynman gauge. The 
renormalization of Wilson loop is performed in the MS— scheme using a general procedure 
tailored to the light-like axial gauge. We find that the renormalized Wilson loop obeys a 
renormalization group equation with the same anomalous dimensions as in covariant gauges. 
Physical implications of our result for investigation of infrared asymptotics of perturbative 
QCD are pointed out. 
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1. Introduction 



There exists a unique possibility in non-Abelian gauge theories to formulate the nontrivial dynamics of gauge 
fields in terms of the gauge invariant collective variables [jy 

Wc = ^(0| trTPexp (ig £dx» Afa)t a ^ |0), 

called the Wilson loop expectation values. Here, A®(x) is the gauge field, T orders gauge field operators in time 
and V orders generators t a of the SU (N) gauge group along a closed integration path C. With this definition, 
Wc is a gauge invariant non local functional of a gauge field depending on the integration path and satisfying 
nontrivial loop equations Q] . 

However, a closer consideration shows [|| [| that Wc is not a well defined quantity for some paths C lying 
in 4— dimensional Minkowski space-time. The problem is that for a closed path C having either cusps, or 
self- intersections ||, or segments lying on the light-cone |Q, Wc, when calculated in perturbation theory, is 
divergent even after renormalization of gauge field and coupling constant. 

Perturbation theory expansion of Wc looks like 



W c = 1 + V(i<?) n t dx^---i> dx^0( Xl > ■■■>x n )±-trG fll ... fln (x 1 ,...,x n ), (1.1) 



where G^...^ (x%, . . . ,x n ) is the n— point Green function and the 9— function orders points x\, . . . ,x n along the 
integration path C. 

Let us consider the simplest closed path in Minkowski space-time shown in fig. 1 with parallel segments lying 
on the light-cone. Under integration over, say, X\ and Xi along the path C, the renormalized Green functions 
Gf_ ll ■■■fj, n (x\, X2, ■ ■ ■ j x n ) may have singularities at coinciding points {x\ = x 2 ) and light-cone singularities ((xi - 
X2) 2 = 0), which give rise to very specific divergences of Wc- 

These divergences were considered for a long time Q as an undesirable byproduct in dealing with the loop 
equations. This opinion radically changed after it was found Q that those specific divergences are of the 
utmost importance in perturbative QCD. It turned out indeed that there is an intimate relation between their 
renormalization and infrared asymptotics of hard processes in perturbative QCD. For instance, the Wilson loop 
expectation values calculated along paths partially lying on the light-cone obey renormalization group equations 
(RG) which coincide with the Altarelli-Parisi-Lipatov and the Brodsky-Lepage evolution equations || |||. The 
"bremsstrahlung" function, well known in QED for a long time, is nothing but the cusp anomalous dimension 
of the Wilson loop j§ §. The same function, called velocity dependent anomalous dimension, was rediscovered 
recently within the heavy quark effective field theory J?]] . That is why the investigation of the renormalization 
properties of Wilson loops has not a pure academic interest. 

In a previous paper || the calculation of the Wilson loop along the light-like path of fig. 1 was performed in 
the second order of perturbation theory in the Feynman gauge. The dimensional regularization was used and 
all divergences were subtracted in the MS— scheme. It was shown that the renormalized Wilson loop, being an 
even dimensionless function of the scalar product (n • n*), of the renormalization point p and of the coupling 
constant, has a form 

W§ = W§(j>,g), p 2 = (n-n> 2 , 

and obeys the following RG equation || || : 

PIT + /^TT lo S W C (P, 9) - -2T cusp (g)(log(p 2 + iQ) + log( V + i0)) - T(g), (1.2) 



dp dg 

n and n* being the two light-like vectors 

rc M = (T,0,0,-T), n* = (i,0,0,L), (n • n) = (n* ■ n*) = 0, (n-n*) = 2LT. 

with L, T > 0. The unusual property here is that the r.h.s. being considered as anomalous dimension depends 
on the renormalization point. It implies that Wc does not renormalize multiplicatively. 

Equation ([F^) contains two gauge invariant anomalous dimensions: T cusp (g) and r(^). The first one is related 
to the asymptotic behaviour of the so-called cusp anomalous dimension and has very interesting properties. In 
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particular, T cusp (g) does not depend on the form of the integration path C and is equal at two- loop order to 

iS 

W.) -** + (*)"*(<*(£- <-> 

with a s = f^. Here, C F = ^± and C A = A are color factors of the gauge group SU(N). The anomalous 
dimension r(g) does depend on C and it was found at two-loop order for the path in fig. 1 to be || 



(1.4) 



The RG equation (1.2) was checked in Feynman gauge and gauge invariance implies that it should be fulfilled in 
any gauge with the same anomalous dimensions T cusp (g) and T(g). This is one of the statements we are going 
to verify. 

In the present paper we calculate the Wilson loop expectation value along path of fig. 1 in the light-like 
axial (LLA) gauge: 

n»Al(x) = (1.5) 

with the gauge fixing vector along one of the sides of the integration path at the second order of perturbation 
theory using dimensional regularization with D = 4 — 2e. 

It is well-known that the free propagator of the gauge field is more singular in axial gauges owing to 
the presence of additional "spurious" singularities. In the LLA gauge we are considering, a prescription, the 
Mandelstam-Leibbrandt (ML) causal prescription , has been proposed to consistently M] define it as 

n (h\ — 1 ( k^n u + kyfi^ \ 

We stress that the +i0 prescription only specifies how the integration contour is to be distorted near the pole. 
In particular it is irrelevant in the expression 

(nk) _ l 



(nk) + iO ■ (n*k) 



always in the sense of the theory of distributions; one can also prove that Mandelstam's proposal does indeed 
coincide with Leibbrandt's one, namely 

1 (n*k) 



(nk) + iO ■ (n*k) (nk)(n*k) + iO ' 

The vector propagator contains the famous (kn) = singularity and, as a consequence, individual Feynman 
diagrams contributing to Wq in LLA gauge may contain additional poles at D — 4. These additional poles 
should be compensated in the sum of all diagrams contributing to Wc, in order to ensure its gauge invariance. 
The mechanism of this compensation is not trivial and very sensitive to the prescription one uses to define the 



propagator (1.6). Our purpose it to show that we shall indeed find full consistency using the ML prescription, 
at least up to the second order in the loop expansion. 

Calculating Wilson loop along the path of fig. 1 in the LLA gauge, we will meet three different kinds of 
divergences: specific light-cone singularities of Wilson loops, singularities which are peculiar of the light-cone 
gauge we are using and conventional divergences of the Yang-Mills theory. Since the first two of them originate 
from the same phase space region of gluon momenta, i.e. gluons propagating along the light-cone, an interplay 
is a priori possible destroying their correct compensation. It would mean that the LLA gauge with the ML 
prescription is sick. 

This doubt was raised in a recent paper pj| . The leading singularities of the Wilson loop of fig. 1 were 
calculated using the LLA gauge with the ML prescription in the second order of perturbation theory and the 
following expression was obtained 

W c ~ exp{Ba 2 s C F C A {D - 4)~ 4 + 0((D - 4)- 3 )}, 

B being a suitable constant, which turned out to be different from the Feynman gauge result: B — 0. 
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In what follows we present a complete two-loop calculation of the same dimensionally regularized Wilson 
loop. As dimensional regularization does not spoil gauge invariance, we shall verify that the expression we 
obtain in the LLA gauge for the unrenormalized (but dimensionally regularized) Wilson loop exactly coincides 
with the analogous quantity evaluated in Feynman gauge up to terms vanishing when D = 4, i.e. not only 
leading and non-leading logarithmic terms but up to finite ones. This is a complete check of gauge invariance 
and thereby of the correctness of the ML prescription, which is crucial to this goal. 

Then, as a second step, we shall renormalize Wilson loop expectation value in the MS— scheme according 
to the theory of renormalization in LLA gauge with ML prescription developed in refs.|L2|. We show that it 
indeed obeys the RG equation (1.2) with the same anomalous dimensions (1.3) and ( |1.4| ) as in covariant gauges. 

The paper is organized as follows. In Sect. 2 the properties of free gluon propagator are described and the 
full set of non vanishing Feynman diagrams contributing to Wc is identified. In Sect. 3 we determine Wc at 
the first order of perturbation theory. Second order Feynman diagrams are evaluated in Sect. 4 and summed 
to get the final expression for the dimensionally regularized Wilson loop, recovering the Feynman gauge result. 
In Sect. 5 the Wilson loop renormalization, which presents non trivial features in LLA gauge, is performed in 
the MS— scheme and a check is made of the anomalous dimensions against the corresponding ones in Feynman 
gauge, thereby supporting the renormalization procedure in LLA gauge developed in refs. |12| . Sect. 6 contains 
concluding remarks. Some technical details of our calculations are presented in the Appendices. 



2. Free gluon propagator 



To calculate the vacuum average Wilson loop we use the definition (1.1). The integration path C is shown in 



fig. 1. It lies on the light-cone and is formed by two light-like vectors n and n* we have already introduced. As 



a consequence, in the LLA gauge (1.5) gluons cannot be attached to any segment of the path C going along the 



gauge fixing vector. We parameterize the two remaining segments of the path C as follows: 

Xl (t)=n% t€[0,l]; x 2 (s) =n + n*s, se[l,0]. (2.1) 

Note that parameter s runs from 1 to in order to provide the correct orientation of the integration path. 

There are a lot of Feynman diagrams contributing to the Wilson loop expectation value at the second order 
of perturbation theory. However, the number of diagrams one has to evaluate can be drastically reduced using 
the non- Abelian exponentiation theorem Jl3[ . According to this theorem 

oo oo 
n— 1 n— 1 

where is given by the contribution of with the maximal non-Abelian color factor to the n— th order 

of perturbation theory which is equal to Cp for n = 1 and C'fCa for n = 2. The two- loop diagrams containing 
color factors Cf and CfCa are shown in figs. 2, 3 and 4 where we omitted symmetric diagrams. It turned out 
that, owing to the properties of the gluon propagator with the ML prescription, many of these diagrams give 
vanishing contributions. 

The peculiar features of the free gluon propagator, when using the ML prescription for the "spurious" 
singularity, have been discussed in fl4|| . In particular it has been shown that the propagator behaves in the 
coordinate representation as a tempered distribution at variance with the expression one would obtain adopting 
the Cauchy principal value prescription. In this last case additional infrared singularities would plague the 
component n*n*G MI/ (a;). 

The explicit expressions for the propagator with the ML prescription in the coordinate representation are 



with the causal scalar distribution 



G a p{x) = 5 a pG{x), a = 1,2, 
1 1 



G{x) = - 



4tt 2 x 2 -iO 



and 



G- a {x) = -£-J* d£G(£,x-,xr), G—(x)=2^-£ d£G(£,x-,xr), G +li (x)=0. 
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We have used here the light-cone variables x_, x+ and xt for the components of vector x. Their definitions and 
some useful identities can be found in Appendix A . 

The origin of this reasonable behavior from a mathematical point of view can be traced back to the causal 
nature of the ML prescription in which the "spurious" pole complies with the Feynman poles position in the 
complex energy plane so that no pinches occur under Wick's rotation of the integration contour. From the 
physical point of view it has been shown jl5| that a "longitudinal" ghost enter the theory. This ghost, while 
decoupling in all physical quantities, is responsible for the mild infrared behavior of the gluon propagator. 
Nevertheless, since the propagator is a generalized function, local limits may not exist. Indeed, using the 
expression for the G (x) component one immediately finds that the function 

G{x) = nlntG^{x) 

is divergent for xt = 0. 

Then one uses the dimensional regularization to define it at Xt = 0, as 

2tt~ d / 2 Y(D /2l (rn*) 2 
G{x +1 x T = 0,*_) ee n> lG^{x +1 x T = 0,*_) = 4 _y \_jll Q)D/ 2 - (2-3) 



Of course, it is the "spurious" (kn) = singularity which is responsible for the appearance of a pole in (2.3) for 

It is now straightforward to realize that one- and two-loop diagrams with all gluons attached to the same 
contour side (see e.g. figs. 2a and 3a,b) give zero. Consider first the contribution of the diagram of fig. 2a: 



{ig) 2 C F ^- D f ds I dtnlnlG^{n k {s-t)). 



Using the definition (2J3) we find that it vanishes as 

nlnlG^{n*{s - t)) ~ {n*n*) 2 - D/2 = (2.4) 

for KeD < 4. The diagrams of fig. 3a and 3c also vanish since the same function enters the corresponding path 
integrals. We stress that the property (2.4) is a peculiar consequence of the ML prescription. To show this, let 
us perform again the evaluation of the diagram of fig. 2a in the momentum representation: 

*-° ( l d« f M f dk + dk - dk T r -ik-LV2(s-t) 2i k- 



( ig) c F » - j Q ds j Q dtj {2n - D e - {2k+k _-e T+iQ) {k+ +mk-y 

It is important to note that fc+ does not enter the exponent when the gluon is attached to the same contour 
side. As a consequence, the integral over the k + component is zero because with the ML prescription poles of 
Feynman and "spurious" denominators lie on the same side of the real k+ axis. The same arguments can be 
applied to show that the diagram of fig. 3b gives zero. In this case we have an integral over two momenta which 
vanishes after one closes the integration contours over "+" components of both momenta without encountering 
any pole. 



3. One-loop calculation 



In the lowest order of perturbation theory, Wilson loop (1.1) has a single nonvanishing contribution coming 
from the diagram of fig. 2b 



i-D 



dtn*nlG^ \x 2 (t) - Xl (s)), 



where the functions xi(s) and X2 (t) are defined in ( |2.l|) . Notice, that the vector xi(s) — x%(t) lies in the 
plane of the Wil son loop and one can use expression ( |2.3| ) for the propagator in Minkowski space-time. After 
substitutions of (2.3) and (2.1) one gets 



W. 



(i) 



r D/2 



,_ D 2T{D/2) i 
4-D 



{nn*f ds dt[2{nn*){s-t)+iO]- D/2 . 



1 In the Feynman gauge, the propagator of the gauge field is equal to — 4 ^T /2 r(Z3/2 — 1)(— x 2 + iO) 1 D ' 2 and 
is regular for D = 4. 



5 



The singularity of the integrand for s = t comes from the propagation of the gluon in fig. 2b along the light-cone 
from point n*s to n + n*t. After an integration which is trivial in the region ReD < 2 and a subsequent analytic 
continuation up to ReD < 4, one finds 



w (i) ff 2 r 4IW2-1) 



(- P 2 + z0) 2 - I5 / 2 + (p 2 + z0) 2 - I? / 2 . (3.1) 



We have obtained the same result in the momentum representation, i.e. integrating first over the contour 
and performing the Fourier transform afterwards. Although calculations are more involved in this way, they 
can be done directly in the strip 3 < ReD < 4. 

We conclude that in the first order of perturbation theory the Wilson loop expectation value contains a 
double pole. This expression coincides exactly with the analogous expression in the Feynman gauge. 



4. Two-loop calculation 

Calculating the Wilson loop in the second order of perturbation theory we use non-Abelian exponentiation and 
restrict ourself only to diagrams containing the non-Abelian color factor CfCa- As was shown in Sect. 2, among 
these diagrams only those shown in fig. 4 give non vanishing contributions whereas diagrams of figs. 3a,b,c do 
vanish. 



4.1. Crossed diagram 

The contribution of the "crossed" diagram of fig. 4a to Wc contains the integration of two propagators over the 
path: 

W crLed = {igfC F {C F - C A /2) f d Sl f 1 ds 2 [ dh f 1 dt 2 G(n + - Sl ))G(n + n*(t 2 - s 2 )). 

Jo Jo Ji Ji 

A 1 (2) 

After substitution of (2.3) W cr ' ossed acquires a double spurious pole before the integration over s— and t— parameters 
which are ordered along the path as shown in fig. 4a 

14 crossed ~ F (C F ~ C A / 2) _ {nil ) 

d Sl f ds 2 f dtx f dt 2 [(2(nn*)(s 1 -t 1 )+i0)(2(nn*)(s 2 -t 2 )+iQ)r D/2 . 

JO JQ Jt x 

As in the previous case, we expect to get additional singularities from the propagation of both gluons along 
the light-cone which correspond to the following values of parameters: s\ = s 2 = t\ = t 2 . Indeed, a careful 
integration for ReD < 2 and therefrom an analytic continuation up to ReD < 4 leads to 

wSL, d = fSpWC - c,/2) 4r ^ D y { [(V + .0)-" + („= + ,o)-«] 

+ 4(l-23i^))[(V +< 0)(^ +< 0)]^}. (4,) 

Again the result has been checked in the momentum representation, where the calculation can be directly 
performed in the strip 3 < ReD < 4. 

(2) (2) 

^crossed ^ as a third pole at D — 4. Thus, the leading singularity of W<, r ' ossed is formed by a double pole 
arising from the "spurious" contributions and only one single light-cone singularity. At the same time an 
analogous crossed diagram in the Feynman gauge gives rise to a fourth-order pole which however cancels in the 
sum of all diagrams. It means that the Feynman gauge is more singular than LLA gauge as far as light-cone 
singularities are concerned. This property is well known in perturbative QCD. As we shall show below, all the 
diagrams of fig. 4 have no fourth pole at all. 



G 



4.2. Self-energy diagram 

ram of fig. 4b 

w s% = {w) 2 C F ^- u J ds J dtGx{n + n*{t - a)), (4.2) 
ion the function Gi(fc) is denned as 

Gi(fc) = ln* M Cr(k)U vp (Jt)G>'\k)nl (4.3) 



The contribution of the diagram of fig. 4b to the Wilson loop reads 

rl i-O 

lo Ji 

where in the momentum representation the function G±(k) is defined as 



Here Gu, v (k) is the free gluon propagator (1.6) and Ii vp {k) is the one- loop gluon self-energy operator in the LLA 
gauge. There are no Faddeev-Popov ghosts and H up (k) gets contribution only from the gluon loop.Q 

The calculation of the one-loop self-energy Tl vp (k) in the LLA gauge is a very cumbersome problem. As 
was shown in []l6|| , the self-energy W vp {k) is transverse and it is decomposed into seven independent tensor 
structures some of which depending on the gauge fixing vectors n and n*. Moreover the contributions to the 
self-energy of the tensors involving n p and n* exhibit non polynomial residues at D — 4, whereas the residue of 
the contribution of the "physical" tensor (g up k 2 — k v kp) gives rise in the LLA gauge to the one-loop (3— function, 
since Z\ = in this gauge as is well known. 

Notice that, after substitution of the self-energy into ( |4.3| ), all tensors do contribute to the function G\{k) 
and the residue at the pole differs from the j3— function. 

The one-loop self-energy was calculated in |16| up to 0((D — 4) 1 )— terms. It is important to recognize that 



these terms, which vanish formally for D — 4, might become divergent after integration over the path in (4.2) 
due to light-cone singularities. This is the reason why one is obliged to first calculate the self-energy for an 
arbitrary D, then to integrate it over the path and only eventually perform a Laurent expansion at D = 4. We 
have performed this calculation in the momentum representation and some details are given in Appendix B . 



We present here only the final expression for the function (4.3) 



4(nn*)\ , J Kf (^fc)K) ) n ZD-2{n*k){nn*) n . (nn*) 

Gl(fc) " + 4 \TfcT " — & — J " D-i jfespj 13 ~ 8 — Ji 

where the integrals I\ ... I4 are listed in Ap pen dix B. We have checked that G\(k), after Laurent expansion 
around D = 4, coincides with the expression (L3) calculated starting from the self-energy U vp (k) found in @. 



By performing the Fourier transform, we obtain the function G\ in the coordinate representation 

Gl{x) ~ -s^° Afl r^TW~ 2 \W~W\ r^/ 2 - 1 ) 

r»(2J/2 - 1) f 6 _ 3D-2 2(D-2) \ | _ 



(4-L>) 2 V D-3 
which is valid for xt = 0. 



Inserting (4.4) into (4.2) and integrating over the path, in a region ReD < 2, we get the following contribution 



from the diagram of fig. 4b to the Wilson loop 

t 

(27T) 



T\3-D/2)T(D-3) , 



(4 - D) A {D - 3) V r(5-Z>) 
r 2 (L>/2-l) / 3L>-2 D 



-T\D/2- i; 



(4-L>) 3 (L>-3) V 4(D-1) D-3 



(4.5) 



which can be continued to the "natural" strip 3 < Re_D < 4. Again the same result can be recovered directly in 
this strip by performing the calculation in the momentum representation. 



(2) 

The expression for W^ e L exhibits the third pole com ing from the double pole of the function ( |4.4| ) and the 



single light-cone singularity of the path integral in (42) for s = t 



2 We do not explicitly consider the contribution of the quark loop, which does not entail any new feature 
with respect to the one in Feynman gauge frjjfl. 
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4.3. Three gluon diagram 

The calculation of the diagram with the three-gluon vertex is the most cumbersome one. At variance with the 
previous diagrams where we have integrated free gluon propagators in the coordinate representation along the 
path, here we have to deal with the additional vertex of a self-gluon interaction. In the coordinate representation 
it implies integration over the intermediate point z. Then the contribution of the diagram of fig. 4c is given by 

W^ glaon = (ig) il -C F C A ^- 2D J\ S1 j\s 2 £ ds 3 

d ° z I 7 ^^e lfcl(;r2(si) - 2:) e lfc2(;i:i(s2) - z) e l ' C3(;i:i(s3) - z) y(fci,/c 2 ,fc3), (4.6) 
J 



where the parameters Si are ordered along the path and gluons are attached to the path at points xi(si), 2:2(52) 

d D k — TT d D kj 
{2tv) d — lli=l,2,3 (27T) 1 



and x 2 {s 3 ) defined in (2.1). Here i A, = n«=i 2 3 ti \i> denotes integration over gluon momentum and 



V(kt,k 2 ,k 3 ) = -iT^ vp (ki,k2,h)G^'(ki)G UI/ '(k2)Gpp>(k 3 )n*^nl,n*^, (4.7) 

r^i/p(/ci, k 2 , k 3 ) being the standard three-gluon vertex expression. Note that integration over the intermediate 

point z ensures momentum conservation k\ + k 2 + k 3 = 0. 

(2) 

Our strategy for calculating W 3 _ gluon is the following. Instead of integrating over z, we first work out the 
expression for V(ki, k 2 , k 3 ) in the momentum representation, then perform independent integrations over the 
three gluon momenta in (4.6) and integrate over z at the very end. 

After some algebra we get from (4.7) an expression for V(k\, k 2 , k 3 ) involving four addenda 

V{k 1 ,k 2 ,k 3 ) = V 1 + V 2 + V 3 + V i , (4.8) 

where the expression for V\ contains only one single "spurious" denominator [kn] = (kn) + iO • (kn*): 

2(nn*) f 7 , +, { (k 3 n*) (k 2 n*)\ „ / 1 1 \ . . (km* 

Vl = TTTTTT ( fc i" ) ~ TT — T + {k 2 n*){k 3 n*) —— - j— )+ {k 2 - k 3 )n* 1 



k\k\k\ \ \ [k 3 n] [k 2 n] J \[k 2 n] [k 3 n]J [km] 

the expression for V 2 has two "spurious" denominators 

(nn*) 2 (km*) ( 1 1\ (nn*) 2 ( (k 3 n*) (k 2 n*) \ (nn*) 2 ( (k 2 n*) (k 3 n*) 



V 2 = 



and V 3 and V4 depend only on two momenta 







K) 2 / 






[km]k\ \ 



[k 3 n]k 2 J [km]k 2 k 2 \ [k 3 n] [k 2 n 



v _ (nn*) 2 (k 3 ~k 2 )n* y _ (nn*Y(k 2 - fciK | (nn*) 2 (k 1 - k 3 )n* 
[k 2 n][k 3 n]k 2 k 2 ' ! [km] [k 2 n]k 2 k 2 [km][k 3 n]k 2 k 2 



It turns out that, after substitution of V4 into (4.6) and integration over momenta k 2 and k 3 , we get a vanishing 
result even before integration over the parameters Si. 

Consider the first term in the expression V4. Its integration over k 3 leads to z — £1(33) and the resulting 
expression, being considered as a function of (k 2 n), has singularities coming from two propagators pj^y and 
k 'i^_ iQ which lie on the same side of real axis in the complex (k 2 n)— plane |] Hence the integral over the (k 2 n) 
component of momentum k 2 vanishes. The same argument applies to the second term in V4 and therefore one 
concludes that V4 does not contribute to (4.8). 

Substituting expressions for V\ , V 2 and V 3 into fl4.6| ) one finds that the result of the integration over momenta 
ki, k 2 and k 3 can be expressed in terms of the following three basic functions: 

tt> f \ f d k jfe x 1 

l(x) = J J^jD e {k 2 + io)[k n y 



3 We have used here the identity k 2 = 2(kn)(kn*) / (nn*) — k\ with kx = (ki, k 2 ). 
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(2) 

Let us first consider the contribution of the expression V 3 into Wf' luon . Since V3 does not depend on ki, 
integration over this momentum leads to a 5— function implying z = x 2 [si) — n + n* si and we eventually get 

^g A C F C A fj, s - 2D (nn*y J dsi J ds 2 J J ds 3 - ^) Fi{-n - n*( Sl - s 2 ))F 1 (-n - n*( Sl ~ s 3 )). 
Using the expression (C.2) for the function Fi(x) with the vector x lying in the plane of the path C given in 

(2) 

Appendix C, we integrate over the parameters Si and obtain the contribution to W 3 J g i uon arising from V 3 
g 4 r^D/2-1) f 8r 2 (3- D/2) u 2 , .j-b/j 



(2^)-° (4-D) 4 [ r(5-D 

1 



L>-3 



[(V+iOj^ + ^ + iO) 4 -^^. (4.10) 



This expression has a fourth order pole in the dimensional regularization parameter and contains two different 
structures as an even function of (nn*). 

Calculation of the integrals coming from Vi and V 2 is more complicated. We shall only sketch here the 
evaluation of one of them taking only two terms from V% as an example: 

2(k 2 n*)(k 3 n*)(nn*) ( 1 1 



k\k\k\ \ [fen] [k 3 n] 



Replacing the vertex V by this expression in (4.6), one gets 

J = ig 4 C F C A ^ 2D (nn*) / d Sl / ds 2 / ds 3 (4.11) 

Jo Jo Jo 

x / d D zF 2 (n + n* Sl - z )——(F 1 {n*s 2 - z)F 2 (n*s 3 - z) - F 2 {n*s 2 - z)F 1 {n k s 3 - z)), 
J os 2 ds 3 



where the functions Fi are defined in (h9) and their explicit expressions can be found in Appendix C. A t 
variance with the pre vious case, the integral over the intermediate point z is not trivial and is given by ( C.3 ). 
After substitution of ( C.3 ) into J we integrate over the parameters Si and eventually get: 

J = --A_C F C A T -^^ [(-2(nn*) + *0) 4 ^ + (2(rm*) + z0) 4 ^] 



1 

2-DI 



/ 2 p/2- 2 l-P 2 /2 2 M(/3 2 ,/3 3 ) 
™ 3 l-/?2 P 2 +P 3 



x / dppt ' P3 , a / / > ( 4 - 12 ) 



where d/3 = dPidp 2 dp 3 S{l -pi -fa- Pa), < Pi < 1 and 

M(P 2 ,P 3 ) = Pt D ~ Pt° + (P2 ~ Ps)(p2 + P 3 ) 3 - D - 

The integral over the parameters Pi for an arbitrary D can be expressed in the form of a convergent series, 
which however does not sum to elementary functions. For our purposes it is enough to expand it in powers of 
(4 — D) with the following result: 

x {(1 - C(3))(4 -£>) + (§- 2C(3) + f£)(4 - Df + 0((4 - Df)} , (4.13) 

where £(z) is the Riemann function. 

The calculation of the remaining term s entering the expressions for V\ and V 2 is analogous. For all of them 
we get the final expressions given in ( |C.4|) , which contain similar integrals over the parameters Pi with the same 
function M(p 2 ,p 3 ) defined in fl4.12|). 
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Summing all expressions (4.10) and (C.4) we obtain the contribution of the diagram of fig. 4c to the Wilson 
loop, we multiply it by factor 2 to take into account the symmetric diagram with two gluons attached to the 
lower contour side and get 



W. 



(2) 



3 — gluon 



C f Ca - [{-p 2 + iO)( P 2 +iO)] 



2-D/2 8r 2 (3 - D/2)T 2 (D/2 - 1) 



(2w) D 

+ [(V + iQ) 4 -° + (p 2 + i0) 4 -°] T(D - 3) 



(4 - £>) 4 r(5 - D) 



{A-DY \ 



1 



(2r(D/2-l)r(3- D/2) 



1 r 2 (L>/2-l) 6-L»r 2 (3-L»/2)\ 
D-2 T{5-D) ) 



D-3 T(D-3) 



t 2 1 2 - 5C(3) 7T 4 5 ' 

(4-L>) 2 12 + A-D 4 + 80 + 4~ C() 



(4.14) 



This complicated expression is valid up to terms vanishing for D = 4. 

The contribution coming from diagrams with three gluon vertex has been computed also in momentum 
representation, namely integrating first over the contour and then over momenta. All integrals in this case can 
be directly performed in the strip 3 < KeD < 4. Details of the calculation are reported in Appendix D. The 
final result of course fully coincides with expression 



14) and provides an independent check. 



We have calculated in this section all two-loop diagrams contributing to the Wilson loop expectation value 
with non-Abelian color factor CfCa- To find two-loop unrenormalized Wilson loop using non-Abelian expo- 



nentiation theorem we sum expressions (4.1), (4.5) and (4.14), skip terms with color factor C F in front and 
identify the resulting expression with the exponent in (2.2). 



Notice that unrenormalized Wilson loop is a finite function of D as long as D ^ 4. Trying to give a physical 
meaning to this quantity, we have to take care of poles at D — 4 and to perform the renormalization procedure, 
which eventually will provide us with the two-loop renormalized Wilson loop for D — 4. Renormalization will 
be performed in the next section. 

At this stage we can already settle our first problem, namely whether the LLA gauge with the ML prescription 
is sick in the calculation of the Wilson loop under consideration. We stress that unrenormalized Wilson loop 
for D 7^ 4 is a well defined gauge invariant quantity. The same unrenormalized Wilson loop was calculated in 
H using Feynman gauge. We take here the opportunity of correcting a misprint in the expression for Wf in 
rcf. j|, where the constant j^tt 2 should be replaced by -^tt 4 - 

Combining two-loop expressions for unrenormalized Wilson loop in the axial and Feynman gauges and 
performing their Laurent expansions near D = 4 we find that Wilson loops are identical in both gauges up to 
terms vanishing for D = 4. It means that LLA gauge with the ML prescription works perfectly, at least up to 
two-loop order. 



5. Wilson loop in the light-like axial gauge 

In the previous section we found the two-loop unrenormalized expressions for all independent diagrams con- 
tributing to the light-like Wilson loop. To give a meaning to these expressions for D = 4 we have to define 
the renormalization procedure which, when applied to unrenormalized Wilson loop, removes all divergences 
including overlapping divergences like the one of the gluon self-energy in fig. 4b. 

Calculating diagrams we have met three different kinds of divergences: specific light-cone divergences of the 
Wilson loops, spurious divergences of the LLA gauge and conventional divergences of the Yang-Mills theory. 
All of them have manifested themselves as poles in the dimensional regularization parameter. 

The problem with the renormalization of the first kind of singularities has to do with to the fact that 
Wilson loop is a non-local functional of the gauge potentials. As was shown in ref.g], had the integration 
path neither cusps nor light-like segments, it would be enough to renormalize Green functions and coupling 
constant to remove all divergences from Wc- But since the integration path under consideration has both kinds 
of peculiarities, one should first renormalize the gauge potentials and coupling constant and then subtract the 
remaining specific light-cone and cusp divergences . 

The general theory of renormalization in LLA gauge has been developed in refs. [ji~2| . In particular it has 
been shown there that the vector potential A^(x) renormalizes at all orders in the loop expansion with two 
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renormalization constants Zj, and Z% 

A,{x) = Zl /2 [A*{ X ) (1 - Z^)n^ R (x)}, (5.1) 

with the non-local gauge dependent quantities O(x) being defined in terms of the covariant derivative D ab = 
S ab d^ + gf abc A^ in the adjoint representation and of the field tensor as 

fl a (x) = [{nVy^inn^n^F^ix). (5.2) 

The operator (nl?) _1 is understood as a series expansion in the coupling constant with boundary conditions 
giving rise to the ML prescription. It is immediate to realize that renormalization does not change the gauge 
condition. 

To our subsequent calculations, only the effective part of the zero-th order expansion in the coupling constant 
ofeq. (5.2) 

n-w - 

(n n) 

will be relevant. 

The coupling constant g is renormalized as Z 3 X ^ 2 g R as Z\ = Z3 in this gauge. As a consequence the Wilson 
loop gets the expression in terms of renormalized potentials 



Wc = i(0| trTPexp (ig R j>^[A R {x) - (1 - Z^)n^ R {x)] 



(5.3) 



The renormalization constants Z3 and Z3 at order g are given by fll 



Z > = 1+ 24^4— D» (5 - 4) 

* = ^i^n- <"> 

The extra-term in the Wilson loop (5.3) gives rise at 0(g i ) only to a contribution coming from the cross-product, 
which turns out to be equivalent to the following correction to the gluon propagator 

6G^(x) = jfLJ^G^ix). (5.6) 

In turn, the divergent part at D = 4 of the gluon propagator due to the self-energy insertion of fig. 4b can easily 
be extracted from [ |l6| . It originates from the following two contributions of the self-energy tensor 

'S 2 ^ [11 1 ,2 t ^ , 4 (nk) 



One can easily check that the contribution from the second tensor is exactly compensated by the correction in 
the expression of the Wilson loop due to the presence of f2. Therefore only the term which is proportional to 
the 13— function is to be subtracted and the resulting path integral coincides (up to factor — -^jCaj^) with 
one-loop expression ( |3.l|) 

(-2(nn*) + iO) 2 - D/2 + (2(nn*) + iO) 2 - D/2 } . (5.8) 



(2) _ ff 4 11 „ n T{D/2- 



m'L = ^^c F c A 



(4-D) 



Once this expression has been subtracted from the sum of contributions to the Wilson loop, turning the coupling 
constant g into g R , not all divergent parts at D — 4 are eliminated. The remaining divergences are just light- 
cone divergences of the Wilson loop, which we also have to subtract taking care of overlapping singularities. 
Here we meet some simplifications in the structure of light-cone singularities in LLA gauge. It is indeed well 
known from perturbative QCD that diagrams like the one of fig. 4, having non-Abelian color factor CfCa and 
either being nonplanar (see fig. 4a) or containing gluon self-interaction, do not exhibit overlapping light-cone (or 
collinear) divergences. It means that, in order to renormalize those diagrams, one can simply subtract all poles 
remaining after renormalization of Green functions. Then, using the non-Abelian exponentiation theorem, one 
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realizes that, in order to renormalizc all light-cone singularities of the Wilson loop W, it is enough to subtract 
poles from the exponent w^ n \ which does not contain any overlapping light-cone singularity. 

After having subtracted those remaining poles in the MS— scheme according to this procedure, we verify 
that the resulting renormalized expression satisfies the RG equation (1.2) with the same anomalous dimensions 
of eqs. (1.3) and (1.4), providing a further test of gauge invariance and, at the same time, of the way in which 
rcnormalization operates in the LLA gauge. 



6. Conclusions 

We end by summarizing what we have achieved in previous sections. We have first evaluated in perturbation 
theory at 0(g 4 ) a light- like Wilson loop for the Yang-Mills theory in LLA gauge with the ML prescription for 
the vector propagator. All singularities, the ones related to the ultraviolet behaviour of the theory, the ones 
peculiar of the contour we have chosen (light-like lines and cusps) and the ones related to the gauge choice, 
have been dimensionally regularized and manifest themselves as poles at D = 4. 

As dimensional regularization does not spoil gauge invariance, we have checked that our result does indeed 
completely coincide with the analogous one obtained in ref.Q for the same contour in Feynman gauge. In this 
way a complete test at 0(g A ) of the correctness of the ML prescription in LLA gauge has been achieved. 

Then we have renormalized our result in the MS— scheme. In so doing we explicitly extend the general 
prescriptions given in refs.[^2| to a non-local operator (the Wilson loop). We again recover in a quite non 
trivial way the Feynman gauge findings; in particular both cusp and anomalous dimensions, which obey the RG 
equation and are relevant to the physics of soft radiation, do indeed respectively coincide with their Feynman 
gauge counterparts, providing a further test in favour of general procedures concerning canonical quantization 
and renormalization in LLA gauge as well as of its usefulness in practical calculations. 
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7. Appendices 

Appendix A. Light-cone variables 

For an arbitrary vector x — (xq, x\,X2,X3) light-cone variables are defined as 

x+ = ~^( x o + %a), x- = ^=(x - x 3 ), x T = (xi,x 2 ) 

The following identities are fulfilled for two 4-dimensional vectors x and y: 

x± — x T , x 2 — 2x+x- — x% : (xy) = x+y- + x-y+ — (xt • yr), d A x = dx+dx-efxr 
For gauge fixing vectors n and n* we have 

71+ = 0, n_ = V2T, n T = 0, n* + = V2L, n*_ =0, n* T = 0. 



Appendix B. Useful integrals in the evaluation of the self- 
energy diagram 

We list here the integrals used in the expression G\ (p) . We follow the conventions given in Appendix A for light- 
cone coordinates. In the defining integrals, ML prescription for the spurious poles will always be understood. 

f d D k 1 
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Up) 



d D k 1 



(B.l) 



[2tt) d k 2 {p + k) 2 [k+] 
i (-p 2 - iO) D / 2 - 2 fr(3-D/2) 



(4tt) d / 2 



(2-D/2) 2 



T 2 (D/2- 1) 

r(D - 3) 



i r(2 - d/2 + s)r(-s) ^ p 2 . 



2vri 



7] — tOO 



D/2 -2 + s 



p 2 + iO 



(B.2) 



(27T) 15 k 2 ( P +k) 2 

i r(2- D/2)r 2 (D/2- 1) 2 



^4 (p) 



(4tt) c / 2 



r(D - 2) 



(-p 2 - l o) D / 2 - 2 



(B.3) 



(27r) B /c 2 (p + /c) 2 [fc+] 
i (-p 2 - iQ) D ' 2 - 2 



( iw )D/2 [p 

(p 2 + »0) 



p_r(2-D/2) 



r 2 (D/2- i) 



2[P+] 



r(i - d/2) 



r(D-2) D/2-1 
T{D/2 - 2)r(L>/2) D/2 -3 



T(D - 2) 



-™ 2 - jn^/ 2 - 1 /■-'»+*«' 



i (-p 2 - iO) 



(47T) 15 / 2 

In Eqs. (B.2) and (B.4), n G (0, 1). 



D/2 - 2_ 
T(3 + s - D/2) 



7] — IQO 



(s + D/2)r(3 + s) sin ns \p 2 + iO 



Pt 



s+2 



(B.4) 



Appendix C. Diagram with three-gluon vertex (coordi- 
nate representation) 

Performing calculation of the diagram with a three-gluon vertex we have introduced three functions defined in 
(4.9). Their explicit expressions are 



Fi{x) 
F 3 (x) 



,1-23/2 



(xn*) 



4vr D / 2 {nn*) J Q 

a-D/2 roo 



daa D ' 2 - 2 ( ^l-o/2 e -ia(2x_ a+ -4/e) 



4tt d / 2 J 



daa D / 2 - 2 e~ tax 



in D/2 



{-x 2 + iO) 



l-D/2 



(C.l) 



i (xn*) 
-7 — ^o{xt) 
tt (nn ) jo 



tt (nn*) x — iO 



for arbitrary D— dimensional vector x in the coordinate space. Although some integrals over a— parameters can 
be done here we remain them undone in order to simplify further integration of the product of these functions 
over x. Expression for F\(x) is greatly simplified as the vector x lies in the plane of the vectors n and n*: 



F 1 (x+,x T = 0,x-) 



I nD/2-l)(xnJ_ { _ x2+ . Q)1 _ D/2 



2lr D/2 4 _ D ( nn *) 



(C.2) 



Using flC.2j ) we have calculated in (4.10) the contribution of V3 — term. Calculation of the contribution due 
to V\ and V2 terms involves integration of products of functions ( |C.l[) over intermediate point z. One of these 
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integrals entering into (4.11) can be calculated using (C.l) as 

J d D zF 2 (n + n* Sl - z) (F x (n* s 2 - z)F 2 (n*s 3 - z) - F 2 {n* s 2 - z)F 1 (n*s 3 ~ z)) 

i r(z>-3) rl 



32TT 1 



A-D 



i — Pi 



(C.3) 



where dj3 = d(3idf3 2 d(3 3 5{l — (3\ — (3 2 — (3$), < /3, < 1 and the notation was used for the function $(/3i, (3 2 , f3 3 ) = 
[2(nn*)(Ei=i s l f3 1 -si)+iO] 3 - D The ^-parameters in flCjt ) have appeared after one has performed the follo wing 
chan ge of the integration variables: ai = A/% (i = 1, 2, 3 ) and integrated over A from to oo. Substituting ( C.3 ) 
into ( 4.11 ) and integrating over Sj— parameters we get ( 4.12 ). 

Different terms ente ring into expressions for V\ and V 2 are grouped into 6 items. Contribution of one of 
them is given by ( 4.13 ). The calculation of the remaining 5 term s is analogous. One first rewrites them as 
integrals of a product of F— functions over z, uses relations (C.l) to express the result of z— integration as 
integral over a— parameters and then integrates over s— parameters taking into account their ordering along 
the path. At the end of this procedure substituting ai — A/3.; we get integrals over f3— parameters analogous to 
(4.12). Summarizing, we can write the contribution of 6 items from V\ and V 2 to the diagram of fig. 4c as 



gA -C F C A [(-2(nn*) + i0) 4 -° + (2(nn*) + i0) 4 -°] ^ J} (h 



6An D r * LV v ' ' v v ' ' J (4 - D) 2 
where I\,I 2 ,I 3 correspond to three items in V± and Ii,!^,!^ to three items in V 2 : 



(C.4) 



h = 



h = 



Ja „2-D/2 n D/2-2 l-P 2 12 2 M(p 2j (3 3 ) 

d(jfl ^ 1 -/32 (32 +P3 

(1 - C(3))(4 -£>) + (§- 2C(3) + f£)(4 - Df + 0((4 - D) 3 ) 
l 

, n nl-D/2 r.D/2-2 oD/2-3,, p 2-D/2, „ v 

apPi P 2 /V (i-/ 3 ! )M(f3 2 ,f3 3 ) 



4r(L>/2- 1) 
(4-D) 2 



r(3-D/2) 



/a = 



^1 -, a D/2~2 

I df30 2 1 - D/2 l3 3 D/2 - 31 ^ 



T(D/2- 1) 
2r(L> - 2) 

M(/3 2 ,/3 3 ) 



1-/3: 



£ - §C(3)(4 — D) + ^(4 - £) 2 + C((4 - £) 3 ) 



d[3Pl- D f3%' 2 -*(3 3 D/2 - 2 M{(3 2 ,f3 3 ) = - 



df3/3 2 1 - D/2 /3^ 2 - 3 f3 3 - 1 M(f3 2 ,/3 3 ) 



AT{D/2- 1) 
(4-L>) 2 



r(3-D/2)- 



T{D/2- 1) 
2r(L> - 2) 



2r(D/2-l)r(3-D/2) 2r 2 (3-L>/2) 



(4-U)2 



(4 - D)r(5 - L>) 



[V(3 - D/2) - i/>(D/2 - 1)] 



/« = 



d/3/3 1 1 - c/2 /3 2 D/2 - 2 /3 3 - 1 M( / 3 2 ,/3 3 ) 



2r(L>/2 - l)r(3 - D/2) T(3 - D/2)T(2 - D/2) 



{A-Df 



r(5 - £>) 



i>(3-D/2)-ip(D/2) 



(D/2 -1) 



where V'( z ) — ^ logr(z)/<iz. Note that l\ was used in (4.12). 

Appendix D. Diagram with three-gluon vertex 
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(momentum representation) 



Among the various diagrams, the evaluation of the three-gluon diagram given in Fig. (4.c) is the most delicate. 
In the main text and in the previous Appendix we have shown how the calculation can be performed in the 
coordinate representation. Since this is the only diagram in which the difference between the computation in the 
momentum and coordinate representation is appreciable, it is instructive to sketch the main lines of calculation 
in the momentum representation. 

We start from Eq. (4.4). The trivial integration over dz reproduces the momentum conservation delta 
distribution. Then, wc can write the thrcc-gluon diagram as 



W 3 - g i uon = 

i=i 

where a = 2ig A C A C F ^- 2D , 

Fi{p,q,r) = 
^2(p,q,r) = 



d u pd u qd u r T t {p, q, r) 



(27T) 



2D 



2 2 2 



2 

R] 

2p_ 

b+] 

2<7- 



Y+\ [q+] 



g( P ,q,r)5(p + q + r) = J2 Wi 



| eye. perm, of (p,g,r) j 







ti--r-} 



2r_ 



[q+\ [r+] 



(D.l) 

(D.2) 
(D.3) 
(D.4) 



and Q(p,q,r) denotes the so called "geometrical factor", encoding all the local geometrical properties of the 
loop contour 



G(p,q,r) 



P- 



2L 



dx 



rO ry 
~ x ~ / dy " / 

J2L J2L 



dz~ 



(e 2lp - L -l) 



P-r- 



-2iLp- _ 2) + 



J2iLr_ 



q-T- 



(D.5) 



In the amplitude factors Ti, all the spurious singularities are defined through ML distributions. G{p,q,r), 
instead, does not contain any pole, and all the singularities in the (/_, p- and r_ variables are fictitious. Since 
the Ti (p, q, r) are antisymmetric under the exchange (q, r) — ► (r, q) , the geometrical factor G in the integrand 
(D.l) can be more conveniently rewritten by performing the change of variables (q, r) — > (r, q) in the second 
term of the curly bracket, eq. (D.5). Thus, a successive integration over dPr leads to 



3—gluon 



d D pd D qe 2l P+ T 



(27T) 



2D 



Fi(p, q) 



( e 2ip-L _ 1 ^ e 2 lq _L _^ 

p 2 _q-p 2 q 2 (p + q) 2 



1 + e 



-2iL(p_+g_) 



»=i 

where the amplitudes T\ (p, q) are now defined by 

1 

+ 



{q-{2p 2 +pq)-p-(q 2 + 2pq)} 

{q-(pq - p 2 ) + p-(pq - q 2 )} 



[p+ + q+] 

^2(p,q) = '^—{p-+2q-) 1 
\P+] 

MP,q) = ¥ ^ ] {p 2 .-ql}-^- ] {q-{q-+^-)}. 



(D.6) 



(D.7) 
(D.8) 
(D.9) 
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To cast T\ in the form (D.7), we repeatedly used the so called "splitting formula" 



1/1 



[9+] [p+ + <?+] [p+] V k+\ [p+ + <i+\ 



(D.10) 



which is an identity, in the sense of the theory of distributions, when the spurious singularities arc defined by 
ML prescription. Eq. (D.6) is particularly convenient for the computation of the three contributions in 
momentum representation. In fact, for each of the terms W^ l \ it is necessary to compute only the term with 
the "1" factor in the last square bracket of eq. (D.6), as the second one, proportional to exp[— 2iL(p_ + <?_)], 
can be obtained from the first by the replacement L — > — L. Following these suggestions and using (for instance) 
Schwinger parameterization for Feynman denominators and ML-distributions, one can get the following results 
for W« and 

(1) _. {ALTf- D e-™ D ' 2 cos(7r£>/2) T(D - 4) 

~ ~ la (2ir) D (D-A) 2 

;; f g r(£)/2-l)r(J/2-2) 1 r 2 (3-£>/2) 



T(D-2) {D/2- l)r(5 - D) 

+4r(2-D/2)r(£)/2- 1) I (D.ll) 

(2) (±LTr- D e-™ D /* cos(7r£>/2) T(D A) 

(2ir) D (D-A) 2 

x | r W2 -wD/2- 2 ) _ 2r(D/2 _ 2)r(3 _ J/2) | (D12) 

It should be stressed that, contrary to what happens in the coordinate representation, the expressions and 
W^ 2 ) have been obtained keeping the dimension D in the natural strip 3 < RcD < 4, and therefore they do not 
entail any analytical continuation. 

The computation of is more delicate and require additional care. Following the same instructions used 
for and W^ 2 \ one can get the following integral form for 

w{3) = ia {ALTf^e-^ 2 cos(^/2) TjD - 4) f 1 2_C 

x I"* ^Vl+l-C] D - b [(ri + l-Cr D/2 -{vQ 2 - D/2 ], (D.13) 
Jo 1 > V 

where ip(Q = [( 4 ~ D — (1 — C) 4 ~ D + 1 — 2(]. The integral in (D.13) presents nasty singularities in ( = 1 if D lies 
in the natural strip. To overcome this problem yet remaining in the 3 < RcD < A region, we can for instance 
perform first an integration by parts in dr\ and write 



Jo 1 "+" V 



D/2-2 



2- D/2 

2 F 1 (2-D/2,3-D/2;A-D/2;C) 
(2- D/2)(3- D/2) 



(BAA) 



f 

Jo 



dv -(v + i-0 D - 5 v 2 - D/2 = r2( f ~^ 2) (i-C) D/2 



1 + ?7 W ' r(6-D) 

x 2j F 1 1 (1,3- j D/2;6-D;C). (D.15) 



Notice that the two hypergeometric series in the r.h.s. are analytic functions in \(\ < 1 if D lies in its natural 
strip. Therefore, the remaining integral in d( can be evaluated integrating term by term in the series expansion 
of the hypergeometrics. The final result is 



(3) . (ALTf- D e-™ D l 2 co^D/2)Y(D-A) f 

(2ir) D (D-A) 2 [D/2-2 



T 2 (3-D/2) 

r(5 - D) 



16 



-3r(3 - D/2)T(D/2 - 1) + 2 



T 2 (D/2- 1) 



r(£>-2) 

tP(n + D/2) - V(n + 3 - D/2) + 



+ 



E 



r(n + 2-D/2) 
r(2 - D/2) ^ (n + 3-D/2)n! 

2 



1 



(ra + D/2- l)(n + £>/2) 
r(n + D/2-l)r(5-£>) 



n + 3 - D/2 



r(4 + n-£)/2) 



(4 - D)T(3 - D/2) 



T{n + 3-D/2) 
- n T(n + 6-D) 



E 



n=0 



r(D/2 - 2) 



T(n + 5-D) T(n + 2) \ 



r(n + 3-D/2) T(n + D/2)J 



+ 2 



r(n+l)r(D/2) r(n + 5-D)r(D/2-l) r(n + l)r(3 - D/2) 



I> + l + D/2) 



T(n + 4-D/2) 



T(n + 4-D/2) 



(D.16) 



where ip(z) = d \ogT(z)/dz. Eq. (D.16) concludes the calculation of and therefore of the whole three-gluon 
graph. Notice that all the series in eq. (D.16) are convergent for 3 < ReD < 4. 

Although not straightforward, it can be checked that collecting all the expressions for IfW and performing 
the Laurent expansion around D = 4, the residues at the D = 4 poles as well as the finite term exactly coincide 
with the corresponding quantities evaluated in the coordinate representation, eq. (4.12). 
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Figure captions: 



Fig. 1: Integration path for Wilson loop. Parallel segments belong on the light-cone. 

Fig. 2: One-loop contributions to the Wilson loop expectation value. The graph (a) vanishes due to causal 
properties of ML prescription. 

Fig. 3: Diagrams containing the following subgraphs give vanishing contributions to Wilson loop expectation 
value. 

Fig. 4: Relevant diagrams contributing to the Wilson loop to the second order: "crossed" diagram (a), 
"self-energy" diagram (b) and "three-gluon" diagram (c). 
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